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Abstract 

We develop a computation model for solving Boolean networks that 
implements wires through quantum ground-state computation and imple- 
ments gates through identities following from angular momentum algebra 
and statistics. The gates are static in the sense that they contribute 
Hamiltonian and hold as constants of the motion; only the wires are dy- 
namic. Just as a spin 1/2 makes an ideal 1-bit memory element, a spin 1 
makes an ideal 3-bit gate. Such gates cost no computation time: relaxing 
the wires alone solves the network. We compare computation time with 
that of an easier Boolean network where all the gate constraints are simply 
removed. This computation model is robust with respect to decoherence 
and yields a generalized quantum speed-up for all NP problems. 

1 Introduction 

The prevailing approach to quantum computation evolved from classical re- 
versible algorithmic computation (Bennett 1979, Fredkin and Toffoli 1982), 
where a stored program drives a sequence of elementary logically reversible 
transformations. In reversible-algorithmic computation a time- varying Hamil- 
tonian drives a sequence of unitary transformations (Benioff 1982, Feynman 
1985). It was then found (first by Deutsch 1985) that entanglement, interference 
and measurement yield in principle dramatic speed-ups over the corresponding 
classical algorithms in solving some problems. 

In spite of this important result, this form of computation faces two possibly 
basic difficulties. Its speed-ups rely on quantum interference, which requires 
computation reversibility. Decoherence may then limit computation size below 
practical interest. Only two speed-ups of practical interest have been found so 
far (factoring and database search), and none since 1996. 

Reversible-algorithmic computation is not the most general form of quantum 
computation. Its limitations justify reconsidering quantum ground-state com- 
putation (Castagnoli 1998, Farhi et al. 2001, Kadowaki 2002, among others), 
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a formerly neglected approach still believed to be mathematically intractable. 
Quantum ground-state computation evolved from classical ground-state compu- 
tation (Kirkpatrick & Selman 1994, among others), a well-developed approach 
competitive with algorithmic computation for solving Boolean networks. A 
Boolean network is a set of nodes (Boolean variables) variously connected by 
gates and wires that impose relations on the variables they connect (Fig. 1). A 
Boolean assignment satisfying all gates and wires is a network solution. Roughly 
speaking, all NP problems can readily be converted to the problem of solving a 
Boolean network. 

In quantum ground-state computation, one sets up a quantum network 
whose energy is minimum when all gates and wires are satisfied. In quan- 
tum annealing, one form of ground-state computation, coupling the network 
with a heat-bath of suitably decreasing temperature relaxes the network to its 
ground state, a mixture of solutions (we assume with no significant restriction 
that there is at least one). Measuring the node variables (Hermitian operators 
with eigenvalues and 1) yields a solution. 

It is believed that quantum annealing yields a (still ill-defined) speed-up 
over its classical counterpart. Quantum tunneling reduces the risk that the 
network, in its way toward the absolute energy minimum, remains trapped in 
local minima (e.g. Kadowaki 2002). However, long simulation times seriously 
limit research on this approach. 

Here we develop a hybrid mode of computation. We implement wires by 
ground-state computation. We implement gates as algebraic identities resulting 
from quantum symmetries and statistics. 

We show that relaxation-computation time is comparable with that an easier 
(more loosely constrained) logical network where all the gate constraints imple- 
mented by quantum symmetries are removed. The comparison is based on a 
special projection method. We show that the relaxation of the actual network 
can be obtained as a special projection of the relaxation of the easier comparison 
network. This projection method shortcuts mathematical complexity and sheds 
light on the nature of this form of computation. 

We conjecture that for this computation mode all hard-to-solve (NP) net- 
works become easy (P) and support this conjecture with plausible estimates. 
Dccohcrence is not expected to be as serious a problem for this computation 
model as for algorithmic computation since the network state is intentionally a 
thermal mixture during most of the computation. 

This discussion of quantum computation still belongs to the realm of princi- 
ples, like other literature on quantum ground-state computation, while algorithmic- 
reversible computation is now almost a technology. Nevertheless it is worth 
starting over with a new approach that might overcome fundamental limita- 
tions of algorithmic computation. 



2 



2 Computation model 



We use a network normal form composed just of wires and triodes (Fig.l). Each 
triode r — properly a partial gate — connects three nodes labeled rx, ry, tz 
(replaced by collective indices in Fig.l) with the sum-2 relation 



q TX + q Ty + q TZ = 2, 



(1) 



where q's are Boolean variables and + denotes arithmetical sum. The three 
solutions are the rows of Table I. 

Each wire w(i,j) — w(j, i) is an equality relation qi = qj between two nodes 
i,j (Table II). 

The example in Fig. 1, with Q = 6 nodes, W = 4 wires (lines), and T = 2 
triodes (dashed triangles), has just one solution: g 3 = q 5 = 0, q\ = q<z = q± = 
<?6 = 1- 
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Fig. 1. A network 



Table III 



In the following we give idealized physical models of the various network 
elements. 



2.1 Nodes and triodes 

Network nodes represent qubits, which here we consider in a most general way 
as commuting Hermitian operators with eigenvalues and 1. In most compu- 
tation, the relations to be satisfied are achieved by a dynamical development. 
In principle we may model any network relation by any valid physical relation, 
however. Here we model the triode relation by a spin identity. In the present 
computation model, each node belongs to a triode r, a spin 1 system, which 
might be two spins 1/2, |cr r i, 5OV2 in units K = 1, in a triplet state, with 
total spin vector s T = ^ (<t t i +0Y2), s TZ = ±1, 0. For each spin 1 we define three 
qubits, 

_ 2 2 _ 2 fC)\ 

qrx — s TX , q Ty — s Ty , q TZ — s TZ1 (Z) 

each representing a node of triode r. 

By the composition of angular momentum, the three qubits of each proton 
pair satisfy the XOR gate equation (Table III). The four rows of Table III 
correspond to the singlet and the three triplet states of proton pair r, spanning 
the Hilbert space TL^ . We use = <S>r=i as tne network space. 
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Note that network nodes are not divided into inputs and outputs of the 
computation process as in algorithmic computation. They are all simultaneously 
present in the network as commuting Hermitian operators, related by time- 
independent identities that we use as gates. 

To simplify the physical model, we assume that the spatial wave function of 
each proton pair r is frozen throughout the computation in a stable antisym- 
metric ground state, like that of the protons of an orthohydrogen molecule. In 
the triplet state the three qubits of any triode r obey the triode sum-2 relation 
of (Q) and Table I. This follows from angular momentum composition and triplet 
symmetry, both extradynamical relations; the triode Hamiltonians are zero in 
n N . 

Let H { r ] be the space spanned by the three triplet states of triode r. T~t N = 
(S)r=i ^ C Tift is the network subspace with all the triodes satisfied. 



2.2 Wires 

We define the frustration Hamiltonian of wire w(i,j) in Tiff by = g (qi — q.j) 2 , 
where g is a coefficient to provide the dimension of energy. The eigenvalues of 
H^j are zero when the wire is satisfied, g when it is not (the wire is then "frus- 
trated"). All commute. Therefore the network frustration Hamiltonian 
is 

H^=gJ2(Qi-Qjf, (3) 

M 

where {w} is the set of all wires in the network. 

Let X T be the exchange operator for the two protons of triode r. H$ is 
symmetric under all the X T (X T iily = H^' X T ), since the q's are. Therefore 
the triplet symmetry projection operator T T of the triode r is a constant of 
motion of H^ 1 . If the initial network state is in nP, under it remains 
in it. The ground state of H^' (in Tiff) hosts a mixture of network solutions, 
since all wires and triodes are satisfied. 



2.3 Ising model 

Hffi is quadrilinear in the spin components s TXi y tZ . As a step toward implemen- 
tation, we show that H^ can be represented by pairwise spin-spin interactions 
by adjoining two "idlers," auxiliary spin-1/2 variables Cj, <jj, to the spin-1 vari- 
ables Si 7 Sj already defined for each wire w(i,j). For convenience we normalize 
the idler variables to eigenvalues Oi = 0,1. One of many suitable frustration 



Hamiltonians for wire w(l,2), connecting the nodes q\ 
collective indices 1, 2) is the bilinear form: 



sf and q2 = a\ (using 



H 



«>(i,2) — 9 



(si + s 2 ) 2 + 5 (si + s 2 ) (-(7i + u 2) + <ti<t 2 + 6(cti + a 2) 



(4) 
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This is chosen so that (as is readily checked) the ground-state projection opera- 
tor of the wire (= nodes + idlers) Po( s i: s 2, &i, 02) nas ener gy eigenvalue and 
includes eigenvectors \s\, s 2 , 01, 02} with all five combinations of spin-1 eigenval- 
ues that satisfy the wire, namely (si = ±1, S2 = ±1) and (si = 0, s 2 = 0), with 
correlated idler eigenvalues si,s 2 . All the states where the wire is frustrated 
(si = ±1, s 2 = 0; si = 0, s 2 = ±1) have energy > g for any values of a\,a 2 . 

The reduced state of the nodes in the ground state is the trace over the 
idlers, 

Po(si, s 2 ) = tr CTl;(T2 p (si, s 2 , 01, ct 2 ) (5) 

Since all the wires are satisfied in the network state Pq(si, s 2 , <ji, a 2 ), they are 
satisfied in the reduced node state p (si,s 2 ). After the network relaxes to the 
ground state, we can find a solution to the network problem by simultaneously 
measuring both node bits qi,q 2 , ignoring the idlers. 

It is then straightforward to construct the network purely out of spins 1/2 
with pairwisc coupling, as in the Ising model. We leave the idlers alone but 
replace each spin 1 by the sum of two spins 1/2 with a coupling that favors the 
triplet (parallel) state over the singlet (antiparallel) overwhelmingly. 



2.4 Heat-bath and coupling 

It is convenient to use heat-bath quanta that are distinguishable from the net- 
work quanta. We use a photon-filled cavity with Hilbert space 1~Lb- ^-b ■ — 
Tiff ® Ti-B is the "system" (=network+bath) space, Ti^ := Tiff ® Hb is the 
subspace with triplet symmetry, all triodes satisfied. 

We denote by Hg\t) the heat-bath Hamiltonian and define the network- 
bath coupling in TL^ by 

#i 4) W = i?E [^(*) • *n + B T (t) ■ 0V2] • (6) 

T 

This couples each proton spin to the small random Gaussian time- varying mag- 
netic field B T {t) of the photon field at the site of the spin. To maintain the 
symmetry that we have assumed, we assume that the two protons of the same 
triode r experience the same magnetic field and that the spatial wave functions 
of different proton pairs do not overlap. 

Therefore triplet symmetry (satisfaction of all triodes) is a constant of motion 
of ffj 4) (t), thus also of the system Hamiltonian H^\t) = + H { ^\t) + 

Hj 4 \t); in fact is already symmetric. 



2.5 Network relaxation process 

With a suitable time- variation of B T (t), H^> relaxes the network to its zero 
point. 

Let \ip,t) be the state of the system at time t. The development of \ip,t) 
is generated by H^ 4 \t) according to the Schrodinger equation. The relax- 
ation of the network state is described by the statistical operator p N (t) := 
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tiB (\ip,t) (ip,t\), where tr^ means trace over the heat-bath degrees of freedom. 
If p N (t) starts in H N it remains in it, since both the wire frustration Hamilto- 
nian and the heat-bath coupling are symmetric. 

A direct estimate of relaxation time is likely mathematically intractable, and 
a simulation is very long. We take a shortcut that also sheds light on the nature 
of this hybrid computation. 

We compare the network relaxation time with that of an easier network 
obtained by replacing all triodes (Table I) by XOR gates (Table III): as if pro- 
ton indistinguishability were suspended - each proton pair were replaced by a 
dcuteron. 

The restriction to Ti^ vanishes: a network of XOR gates and wires is loosely 
constrained and easy to solve. In particular qi = for all i is always a solution. 
A XOR network is solvable in poly(Q) time in classical computation and, one 
reasonably supposes, in the present hybrid computation also. Note that the 
XOR gates of the comparison network are also extra-dynamical; they represent 
a physical law, namely the composition of angular momentum. 

That the relaxation of the comparison network is quick can be plausibly seen 
as follows. 

We first replace Hffl by the new Hamiltonian 



rr'(4) _ „(4) 
H N ~ N 



(7) 



Since each triode has exactly two nodes equal to 1 (Table I), we have H'^p = 
+ ^^g) ' ^ n ^ e case °f the actual network we have merely multiplied 
Hff by a constant factor and the replacement is inessential. The ground states 
of Hffl and H'^p are the same, and the respective energy landscapes are pro- 
portional. 

Not so for the comparison network, no longer restricted to TC N . If g g, 
the energy landscape of H'jp has a gradient everywhere toward the solution = 
allowed by Table III. There are no local minima that can trap the comparison 
network on its way toward the absolute minimum; thus the relaxation time of 
the comparison network is reasonably poly(Q). 

We conjecture that introducing H'^ is unnecessary. The XOR gates, being 
extra-dynamical, would not affect relaxation time; they could be removed, which 
would leave us with a set of independently relaxing wires, and no local minima. 

2.6 Comparison system 

The asymmetric Hamiltonian of the comparison system in is H A (t) = 

H%>+Hg\t)+Hf(t), where 

Hf(t) =9^2 [BtiV) ■ + B T2 (t) ■ a T2 \ (8) 
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is the asymmetric coupling. Now we have two independent random Gaussian 
time- varying magnetic fields at each proton site, such that 



B T {t) = [B Tl {t) + B T2 {t)\ /2 (9) 

is the actual heat-bath. This is always possible since the sum of two Gaussian 
distributions is also Gaussian. 

Let \ip, t) be the state of the comparison system, whose development is gen- 
erated by H A (t). The comparison network relaxation is described by the sta- 
tistical operator p^(t) := tr# (\<p,t) (<p,t\). 

2.7 Continuous projection method 

The symmetrization operator for all the proton pairs is 

T 



It projects on H {3 1 Clearly PH A {t)P = H?\t), and so 

PH A (t)P = H^(t), (11) 

given that H^hs symmetric and P is the identity in TLb- 

The development of the actual system (hard triode network and bath) in 
7^(3) j s driven by _ff( 4 )(f); that of the comparison system (easy XOR network 
and bath) in TL^ is driven by H A (t). We show that, given (11), the continuous 
projection on of the development of the comparison system yields the 

development of the actual system. 

Let \tp, t) be an initial state of the actual system in Tt^ , therefore P\ip,t) — 
\i>,t) . Under H^(t), it develops into 

| V, t + dt) = (l - iH {i) (t)dt) | V, t) . (12) 

Under H A {t) 1 it develops into 

\<p, t + dt) := (1 - iH A (t) dt) |V, t) , (13) 

in general non-symmetric. We restore particle indistinguishability by projecting 
\<p, t + dt) on symmetrizing it: 

P \tp, t + dt) = (P 2 - iPH A (t)P) |V, t) = (l - iH (4) {t)dtj |V, t) = | V, t + dt) . 

(14) 

We can see that the continuous projection of the comparison development yields 
the actual development. 
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2.8 Comparing computation times 

Computation time is by assumption poly(Q) for the comparison easy XOR 
network. To estimate that of the actual hard triode network, we decompose a 
AT into N = AT/ At consecutive time slices Atj = [tj, t i+ {\ of equal length At. 
Within each Ati, we consider the relaxation of the comparison XOR network in 
Ti.( 4 \ described by p%{t). At the end of each Ati, we project p%(t) on TL^\ 
then take the limit At — > 0. This yields the actual network relaxation p N (t). 
Within each Ati we consider the decomposition 

p%(t):=p (t)+p F (t)+p v (t). (15) 

• p (t) describes networks with satisfied triodes and wires, namely solutions 
of the actual network; its probability is po{t) := trp (i). 

• p F (t) describes networks with satisfied triodes and at least one frustrated 
wire; Pf(*) := trp F (i). 

• py(t) describes networks with at least one violated triode, wires are either 
satisfied or frustrated; pv(t) ■= trpy (t) . 

We have considered all the possible states of the comparison network. There- 
fore Po(t) + pir(t)+ Pv(t) — 1- Pv(t) goes to zero with At and is annihilated by 
each projection. 

The actual network-bath interaction soon randomly generates a Po(th), a 
mixture of solutions of the actual network, with an extremely small probability 
Po(th) = O (l/2 Q ). For a given confidence level, th does not depend on Q. 

For t > th we apply the projection method. Po(th) = O (l/2® ) becomes 
the nucleus of condensation of the network solutions. 

Within each and every Ati, we take a constant- average logarithmic rate of 
decrease k of the frustration energy of the comparison network: 

E N (t i+1 ) = {l-kAt)E N (U). (16) 

We will show later that there is no error in taking a constant-average rate. The 
relaxation time constant 1/fc is by assumption poly(Q). 

We have E^(t) :~trp^(t)H^ = tip F (t)H^ . In fact there is no contribu- 
tion from p (f), which is the ground state of , and a possible contribution 
from p v {t) would anyhow be second order infinitesimal. Let Pp\t) be the j-th 

(population) element of the diagonal of p F (t). Of course ^2Pp\t) = trp F (t) = 

j 

PF(t). is diagonal, thus we have i?jv(t) = Pf\^)^^' where is the 

3 

j-th diagonal element of H^' . Therefore E N (t) and Pir(t) go to zero together. 
Thus on average: 

PF {t l+1 ) = {l-kAt)p F {U). (17) 
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The decrease of pF(t) implies an equal increase of po(t) + Pv{t)- It is reason- 
able and conservative to consider the increase of Pv(t) dominant. In fact the 
relaxation of the comparison network is quicker because triodes can be violated. 

Note that we compare relaxation rates, not directions: the comparison net- 
work can head toward Tiff ~ Tiff, the actual network remains in Tiff . It is 
like comparing the speed of the keel and the wind on a broad reach. Speeds are 
proportional, while keel and wind go to different places. 

Furthermore, Hff does not couple p (t) with p F (t) or py{t). In fact Hffp Q (t) 
= p Q (t)Hff = 0. Therefore po{t) neither decreases nor increases on average. 
Since PF^t) decreases andpo(i) does not, the ratio PF(t)/po(t) decreases. When 
we project on Tiff at the end of Atj, we remain with a smaller PF(t) and a 
larger po{t) (probability of solutions of the actual network). 

We can focus on the "take-off' of the probability of solution from the ex- 
tremely small value po(th) = O (l/2 N ) to Po(t) close to 1, say po{t) — 1/10. 

During take-off and within each AU we have Pf(£) ~ 1; Pv(t) grows from 
Pv(U) = to pv(ti + At) = kAtpp(ti) « kAt, because of (17) and the assump- 
tion that po(t) remains unaltered. The projection at the end of AU annihilates 
Pv(U + At) reducing p^(t) by about kAt. Renormalizing p%{t) then multiplies 
Po(t) by about (1 - kAty 1 w 1 + kAt at each At*. After a time AT = TV At 
and in the limit At — > 0, we obtain for the actual network: 

Po (t h + AT) » Po (t h ) A lim o (1 + kAt)^ = Po (t h )e kAT « ^e kAT . (18) 

The probability of having solutions of the actual network becomes O (1) in a 
time AT w Q/k = Qpoly(Q)=poly(Q). 

Using a different fcj for each At i; with average value k = ^ i k i At i / 'AT, 
yields the same result: e kAT in (18) should be replaced by fJi e fe,Atl = e kAT . 

3 Conclusions 

The extradynamical algebraic relations expressing particle statistics and angular 
momentum composition can replace the dynamical algebraic relations following 
from equations of motion as computational gates. In this new form of quantum 
computation, the gates of a Boolean network are always satisfied as constants of 
the motion, leaving only equality relations (wires) to be implemented dynami- 
cally. This form of quantum computation is expected to be robust, since it relies 
on thermal mixtures, not pure states, and is plausibly conjectured to be fast, 
turning all NP problems in principle into P. As in quantum algorithmic compu- 
tation, the speed-up is due to the extradynamical character of the computation 
(Castagnoli & Finkelstein 2001, 2002). 

This model of computation highlights the conceptual difference between how 
structures can be assembled in the classical and quantum domain. Quantally 
it is as though one could assemble a jigsaw puzzle simply by piling the pieces 
up and letting gravity lower them into mutual positions that solve the puzzle, 
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analogously to quantum wire relaxation. Classically this way of assembling the 
pieces would be plagued by local energy minima. One may wonder whether the 
assembly of biological molecules under hydrophobic pressure draws on similar 
quantum effects. 

Many of the ideas propounded in this work were developed through discus- 
sions with Artur Ekert. 
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